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A NOTE ON THE ENUMERATION OF DIRECTED ANIMALS VIA 

GAS CONSIDERATIONS 
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In the literature, most of the resuhs about the enumeration of 
directed animals on lattices via gas considerations are obtained by 
a formal passage to the limit of enumeration of directed animals on 
cyclical versions of the lattice. 

Here we provide a new point of view on this phenomenon. Using 
the gas construction given in [Electron. J. Combin. (2007) 14 R71], 
we describe the gas process on the cyclical versions of the lattices as a 
cyclical Markov chain (roughly speaking, Markov chains conditioned 
to come back to their starting point). Then we introduce a notion of 
convergence of graphs, such that if (Gn) — > G then the gas process 
built on Gn converges in distribution to the gas process on G. That 
gives a general tool to show that gas processes related to animals 
enumeration are often Markovian on lines extracted from lattices. 

We provide examples and computations of new generating func- 
tions for directed animals with various sources on the triangular lat- 
tice, on the %, lattices introduced in [Ann. Comb. 4 (2000) 269-284] 
and on a generalization of the lattices introduced in [J. Phys. A 
29 (1996) 3357-3365]. 



1. Introduction. Let G = (y,E) be a directed graph with set of vertices 
V and set of oriented edges E. Let A and S be two subsets of V, with S C A. 
We say that A is a directed animal (DA) with source S if and only if any 
vertex of A can be reached from an element of S through a directed path 
having all its vertices in A (see Figure 1). The vertices of A are called cells 
and the number of cells, denoted l^j, is the area of A. We denote Gs t^ie 
generating function (GF) for DA on G with source S counted according to 
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Fig. 1. Example of a DA with area 6. The cells of the DA are dark and the vertices of 
the source double circled. 

their area 



where is the number of DA on G with source S and area k. 

In the following, we will always assume that the cells of S form an in- 
dependent set on the directed graph G — we say that is a free set — the 
formal definition follows. 

Definition 1. Let G = {V, E) be an oriented graph and x and y be two 
vertices of G. We say that x is a father of y or equivalently that y \s a. child 
of X if there is an edge from x to y. 

More generally, x is called an ancestor of y if there exists a directed path 
from X to y. 

Let now S" be a subset of V; we say that S" is a free set of vertices of G if 
and only if for every x,y £ S such that x ^ y, x is not an ancestor of y. 

In this article we focus on the link between enumeration of DA and hard 
particle gas models. 

Definition 2. Let G= {V,E) be a graph, a gas occupation or gas con- 
figuration on G is a map X from V to {0, 1}. The vertices v £V such that 
X{v) = 1 are said to be occupied, the others are said to be empty. 

A hard particle gas occupation of a graph is a gas occupation with the 
additional constraint that two occupied vertices cannot be neighbors (the 
occupied cells form then an independent set). 

A gas model is a probability law on gas occupations. For a given gas 
model, we call density in a vertex v the probability for v to be occupied. 



A, DA fe>|5| 



with source S 



that is F{X{v) = 1). 



Since the pioneering work of Dhar [7], the connection between DA and 
gas models have been widely exploited. We shall now give a short overview 
of the different contributions on this subject (we refer the reader to [4] 
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Fig. 2. The same DA on two representations of the cylindric version of the square lattice 
with a width of 6. 

and [10] for more exhaustive references). In Dhar [7], using some statistical 
mechanics shows that computing the area generating function for DA on the 
square lattice is equivalent to computing the density of a hard particle gas 
model. This result was obtained after Nadal, Derrida and Vannimenus [11] 
and Hakim and Nadal [9] obtained the generating function of DA on some 
"cylindric" square lattices. 

Those "cylindric" lattices are defined as follows: Let G be an oriented 
lattice — that is, an oriented translation- invariant graph — with its vertices 
indexed by a subset of 1?. If we consider that the abscissa of the vertices 
of G are labeled by elements of [N] := TLjNTL instead of Z, we obtain the 
width-bounded variant of G with cyclic boundary conditions (see Figure 2). 
We denote it G^^^ and call it the cyclic or cylindric version of G of width 
N. 

Bousquet-Melou [4] extends Dhar's correspondence between the hard par- 
ticle gas models and enumeration of DA on cyclic square lattices. Particu- 
larly, she shows that gas models allow the enumeration of DA not only 
according to their area but also for instance according to their left perime- 
ter or their number of loops. Those results were then generalized to a family 
of lattices in a joint work with Conway [5]. In Conway [5] and Bousquet- 
Melou [4], the gas models studied are defined on the cylindric versions of 
graphs and the GF for DA is obtained as the formal limit of the density of 
the gas when the width grows to infinity. Since computing the density of 
the gas model is not always tractable, the former result does not necessarily 
lead to effective results about enumeration of DA. However, that new link 
establishes gas models as a powerful and polyvalent tool for the counting of 
DA. 

In the latter works, the link between DA and gas is formal and appears 
because DA and gas models are shown to verify the same recursive decom- 
position along with the layers of the graph. It notably implies that that 
approach is only valid for graphs that can be decomposed nicely into layers. 
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Le Borgne and Marckert [10] give a new insight into the connection be- 
tween gas and DA. They construct a couphng between random DA and 
random gas models and give a combinatorial proof that for a free set S the 
GF of DA with source S is equal to the probability for the vertices of S to 
be occupied (a construction of that coupling is sketched in Section 2). Con- 
trary to the construction on cylinders, in [10] the gas model is well defined 
on any acyclic graph and in particular on the whole lattice, where some of 
its stochastic properties can be studied. On the square lattice for instance, 
its restriction to a line is shown to be Markovian, which allows to compute 
explicitly the GF for DA with any source included in a line. 

We must mention that there exist other fruitful approaches to the combi- 
natorics of DA. Some results have been obtained by establishing links with 
heaps of pieces introduced by Viennot in [12] (see for instance [2, 3, 6, 13]) 
or with paths in the plane [8] or via the ECO method [1]. 

We now describe the content of this paper and its organization. Our aim 
here is to give a general framework that allows to reduce the enumeration of 
DA with various sources on a graph G to the same enumeration on "simpler" 
versions of G. As mentioned above, simplifying the graphs we work on is a 
classical idea. Here, the difference with the works cited above relies on the 
fact that thanks to the gas construction given in [10], we can now study the 
convergence of the gas models as stochastic processes and not only the formal 
convergence of their density. This leads both to a better understanding of 
the gas models and to new results about enumeration of DA with various 
sources. 

The first point is to make the notion of "simpler" versions of G accu- 
rate; in Section 3.1 we provide a distance on the set of graphs with marked 
vertices, corresponding to sources [see equation (3.1)]. Roughly speaking, 
Gn converges to G for the notion of convergence of graphs induced by that 
distance (which corresponds roughly to the convergence of the neighbor- 
hood of sources) implies that converges to G^. In terms of probability, 
that means the convergence of the finite-dimensional laws of the gas under 
additional assumptions (Theorem 2). 

Then we need to compute the law of the limiting gas process obtained 
thanks to that convergence. That is possible on some lattices. The multi- 
plicative formula obtained for the distribution of gas restricted to a line 
on the cylinder in [4] and [5] leads to the intuition that that multiplicative 
structure may be preserved when the width of the cylinder goes to infinity 
and that the limiting process obtained above should be Markovian. For that 
reason, in Section 3.2 we define a cyclic Markov chain as a Markov chain 
conditioned to come back to its initial state after a fixed number of steps 
(Definition 6). We then give a representation of the gas on the cylinder as 
a cyclic Markov chain. Then in broad terms when the width of the cylin- 
der grows, the conditioning induced by the cyclic condition is less and less 
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constraining. At the limit, it eventually disappears which therefore yields 
that the limiting process is Markovian. We provide in Theorem 3 a formal 
statement of those two ideas; that provides a frame in which the gas process 
on a line is Markovian. 

We define in Section 4.1 the family of lattices (>Cij)/jcN) which extends 
the family of lattices {Cn)n>2 introduced in [5]. We apply Theorem 2 and 
Theorem 3 to it, to the triangular lattice and to the family of 7^ lattices 
introduced in [6]. In Section 4, we show that for those three examples, the 
restriction of the gas process to a line is Markovian. Thanks to the link 
between gas models and GF of DA, that allows us to obtain some GF for 
DA with various sources; see, for example. Proposition 1 for some results on 
the triangular lattice. 

2. Definition of the gas model. We sketch the construction of the gas 
model given in [10] and its link with enumeration of DA according to their 
area. Let G = iV, E) be a directed graph without multiple edges nor directed 
cycles and such that the number of children of each node is finite. 

The probability space we work on is 17 = {a, 6}^ endowed with the a- 
field generated by the finite subsets of vertices. We equip that space with 
the product probability Pp = {p5a + (1 — p)5h)®^ : where 5a is the standard 
Dirac measure on {a}. In other terms, w G is a coloring of G and under 
Pp each vertex has, independently of the others, color a or h with respective 
probabilities p and 1 — p. For x £V, u>{x) gives the color of x. From that 
random coloring we construct DA and a model of gas. Notice that the DA 
and gas process defined below are deterministic functions of the random 
coloring. 

Definition 3. Let S be a subset of V and w be a random coloring of 
G. We denote by S,{io) = {x £ S,uj{x) = a}, the (random) subset of S with 
color a. We then define the random variable A'^ as the maximal DA for 
the inclusion partial order with source S,{uj) and set of cells the a-colored 
vertices x that can be reached from S,{uj) by an a-colored path (see Figure 
3). 

For a set S such that \S\ > 1, the random DA A"^ may be infinite with 
positive probability. Let p^j^ be the threshold for the existence of an infinite 
DA with positive probability (it corresponds to the critical probability for 
the oriented percolation on G) 

(2.1) := sup{p, S : Pp(| A"^| < cxd) = 1 and \S\ < oo}. 

For a general graph G, p^j^- is difficult to compute and can even be equal 
to zero. In the examples given in Section 4, the outdegree of any vertex is 
bounded and in that case pcrit > (see, e.g., Proposition 2.2 of [10]). For 
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coloring of the vertices. Cells colored with a (respectively h) are dark (respectively white) 
and the vertices of S are double circled. 



any p < p^j^ a gas occupation X on G is defined from a random coloring 
u> = {u>{v))i,^v as follows (see Figure 3 for an example): 

rO, '\ioj{v) = h, 

(2.2) X^(^) = | ]-[ [l-X'^iv')), \iu{v) = a. 

\ ti' children of v 

The definition of p'^rn ensures that the gas process is almost surely well 
defined as its recursive computation ends within a finite number of steps for 
any p < p^j,. (see Proposition 2.4 of [10] for details). Prom now on we always 
assume that p < p^^^ and that the gas model considered is the probability 
law denoted induced by that construction. 

The link between enumeration of DA and that gas model is given by the 
following result: 

Theorem 1 (Le Borgne and Marckert [10]). Let G = {V, E) he a directed 
graph and S he a free set of G. For any p in [0, ii^), we have 

(2.3) F^(X«(^;) = l,veS) = (-l)l«lgf (-p), 
where is the radius of convergence of . 

With that theorem, the computation of the generating function for DA 
comes down to the computation of the probability for some vertices to be 
occupied for the gas model Fp. That explains why in the next section we 
focus only on the study of the gas model and resume the enumeration of DA 
in Section 4. 

3. Convergence of graphs, gas models and DA. We develop in that sec- 
tion some tools allowing us to reduce to simpler graphs the study of the 
stochastic properties of a gas model on a graph. 
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3.1. Convergence of graphs. As recalled in the Introduction, most of the 
results obtained about the enumeration of DA on a lattice G via the study 
of gas models have been proved by a passage to the limit. More precisely, 
the gas models are studied on G^^\ the cylindric version of G (see the 
Introduction). For a fixed size n, the set of DA with size n coincide on G 
and G^^'^ when say N n. It amounts to saying that 



or equivalently in the gas model's point of view that the density of the gas 
converges formally (in the sense that J2 o-n,kX^ — >n o-k^^ if and only if 
an,k — >n ak for every n G N). 

The aim of this section is to make clear a notion of convergence of graphs 
(i.e., a topology on the set of graphs) which induces the convergence of the 
finite-dimensional laws of the gas process and hence the convergence of the 
generating function of DA. That convergence is no longer seen only as a 
formal convergence of generating functions but as the convergence of the 
distribution of a stochastic process. 

In the following, we always assume that the graphs considered are di- 
rected, without directed cycles nor multiple edges and that the number of 
children of each node is finite (a node can, though, have an infinite number 
of parents) so that the gas model given in Section 2 is defined. 

Definition 4. We cah marked directed graph, a pair {G = {V,E),Z) 
where Z is a subset of V. We denote by Vz the subset of V of nodes having 
at least one ancestor in Z, and by G{Z) the subgraph of G having as set of 
nodes Vz (and set of edges the edges of E linking them). 

To see Z as a source and G{Z) as the maximal DA on G with source Z 
may help to understand better Theorem 2. 

Definition 5. Two directed marked graphs (G = {V,E),Z) and {G' = 
{V' , E'),Z') are said to be isomorphic — we write (G, Z) ~ [G' , Z') — if G{Z) 
and G'{Z') are equal up to a relabeling of the vertices, in other words, if 
there exists a bijective application (j) from Vz onto V^, such that for any x,y 
in Vz, {x,y) ^ E \s equivalent to (</>(x), 4>{y)) G E' . 

The relation ~ is an equivalence relation on the set of marked directed 
graphs. We denote by O the set of directed graph quotiented by that relation. 
For any marked graph {G,Z) we denote by {G,Z) its class in O. 

We denote A'^ the set of DA on G with source Z. The graph {G, Z) is 
the right (or minimal) structure that provides all the knowledge necessary 
to study the gas configuration on Z and the DA with source Z in G [that 
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depends also on the coloring on {G,Z)]. From the construction of the gas 
model and random DA given in Section 2, it is clear that if {G, Z) ~ {G' , Z') 
then \Z\ = \Z'\ and = Q^i and the application (f) provides a probability 
isomorphism between the gas occupations on Z and Z' , which implies that 
pG(Xf = 1, , G Z) = Pf (Xf = l,s€Z'). 

For any r > 0, we define Br{G,Z) as the subgraph of {G,Z) containing 
only the vertices v of {G,Z) such that d{v,Z) = iniu^z d{u,v) < r, where 
the distance must be understood as a directed distance on graphs, that is, 

d{u,v) = mf{\w\, where w is an oriented path from u to v}. 

As announced above, we now define a distance do on O which gives a 
suitable notion of convergence of graphs: For any O and O' in O, we set 

(3.1) do{0,0')=mil^-^,r such that 5,.(G, Z) ~ B,.(G', Z')|, 

where {G,Z) G O and {G',Z') G O' [we let the reader check that that is 
indeed a distance in O and, in particular, that it does not depend on the 
choices of (G,Z) and {G',Z')]. 

Theorem 2. Let (Gn = (Vn, En), Zn) be a sequence of directed marked 
graphs, and [G = {V,E),Z) be a directed marked graph. Let an,k = G 
1^1 = k} be the number of DA with source Zn in Gn having k cells, and 
denote by Ok = #{A G A^, \A\ = k}. 

If doiiC;^), iG;Z)) ^ then: 

1- Qz!^{p) = Efc>|z„| (^n,kP^ — >n^oo Qziv) = T.k>\z\ "'kP^ ^hcre the conver- 
gence holds formally in the set of formal series with coefficient m N (i.e., 
for any k, On.k cik when oo). 

2. // there exists c,d>0 such that for any n large enough, 

(3.2) On.k ^ cd^ for any k>l, 

then for any p <l/d, the finite- dimensional laws of the gas occupation on 
Zn according to P^" converge towards those on Z distributed according 
to Pp , that is, 

P^" (xf " = i,seZn)^ P^(xf = i,sez). 



Proof. 1. First, if do{{Gn, Zn), {G, Z)) 0, then for any r, when n is 
large enough, the two graphs Br{Gn, Zn) and Br{G, Z) are isomorphic. That 
implies that the coefficients of Q'^'^ and coincide at least up to the rth. 

2. First, condition (3.2) implies that p^-'^. >l/d, therefore the gas model 
P^" is well defined for any p < 1/d. 



GAS AND ENUMERATION OF DIRECTED ANIMALS 



9 



From the construction of the gas model, we can notice that the event 
{X^ = 1, s ^ Z} does not depend on the coloring of all the vertices of G but 
only on vertices of (see Definition 3). Since we assume p < 1/d, is 
almost surely finite according to P^; that implies that for any e > 0, there 
exists iris such that P^(|A^| > mg) < e. 

As when n is large enough, the two graphs Bm^{Gn, Z^) and Bm^{G,Z) 
are isomorphic, there exists an application cj) that maps Bm^{Gn, Zn) onto 
Bm^{G,Z). Thus (j) induces a probability isomorphism between the coloring 
of Bm^{Gn, Zn) and of Bm^{G, Z). Therefore, conditionally on the event 
{|A^| < mg}, the image of A^" by (j) is A^ and we get 

P^"(Xf" = l,s G Zn\\A^\ < m^) =P^(Xf = l,s e Z\\A^\ < m^). 

This concludes the proof, since P^(| A^| < rrig) > 1 — e by definition of m^. 
□ 

Remark 1 . Even if in the applications of the latter theorem in Section 
4, we always assume that Zn and Z are free sets. There is no such assumption 
in the theorem and Zn and Z can be any sets. 

3.2. A variation on Markov processes. The spirit of this section is guided 
by the results obtained for enumeration of DA in [4] and [5] . It often happens 
that the probability distribution of the gas has a multiplicative form on 
cylinders. That leads to the intuition that the limiting process obtained 
when the width goes to infinity is Markovian. We give here an appropriate 
frame to make that intuition rigorous. 

In this section, we always assume that E is a finite state space, v a 
probability measure on E and M a stochastic matrix on E. We say that 
Y = {Yi)i^fq is a (i/, M)-MC if it is a Markov chain with as initial law and 
M as transition matrix. 

Definition 6. For any nonnegative N, we call cyclic Markov chain 
of length N on E with initial law and transition matrix M, a process 
(Xj)jg|Q jv-i} which is a Markov chain conditioned to come back to its 
starting point after steps and we say that X is a (i^, M, A^)-cyclic MC. 

Let y be a (i^, M)-MC, for any xq, . . . ,xn-i £ E, the law of (Xj)jg|o^...^Ar_i} 
is equal to 

P(Xo = Xo, . . . , Xn-1 = XN-l) 

(3.3) 

= FiYo = xo,...,YN-i = XN-i\Yo = Yn). 

In other words. 
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where xn = xq and Zn = Ex[„...,x'^_, ^i^'o) llf=o^ ^xr,x'.^^- 

Note that if X is a (z^, M, A^)-cychc MC, the distribution of Xq is given 

by 

(3.5) P(Xo = x) = for anyxGE 
and the distribution of Xi by 

(3.6) p(Xi = xi) = {y^ T^{xo)m^o,-,{m^-%,^,)i (zn)-\ 

^ XQ ' 

Equation (3.5) imphes that the distribution of is not v except for excep- 
tional cases. Combining equations (3.5) and (3.6) imphes that if v = Ue, the 
uniform law on then the cyclic MC is stationary, that is, for any x € E, 
F{Xi=x)=¥iXo = x). 

On the other hand, assume that the initial law v is an invariant law for M, 
then a (j^, M, A^)-cyclic MC is not necessarily stationary. Roughly speaking, 
the term (M.^~^)xi,xo which appears in (3.6) prevents that probability from 
simplifying even if v is an invariant measure associated with M. 

We now give the main result about the convergence of cyclic Markov 
chains. 

Theorem 3. Let E be a finite state space and \ he a square nonnegative 
matrix indexed by the elements of E such that V admits a simple real eigen- 
value A greatest in modulus than every other eigenvalues. Let {X^^^)n>i be 
a family of stochastic processes such that for every N , X^^^ is indexed by 
{0,...,A^- 1} and 

(3.7) P(xi^)=xo,...,x{f_\ = x,v-i) 



trace(V^) ' 

with the convention xn = xq. 

Let R = {Ri)i(zE and L= [Liji^E be respectively a right and a left eigen- 
vector associated with A such that their dot product is equal to one, that is, 
J2 LiRi = 1. 

(i) For each N>1, X(^) is a {Ue,M, N) -cyclic MC, where M is equal 

to 

(3.8) M,,,- = Vi,,-^^ for i,j G E. 

(ii) Let now X = (Xj)jgN be a (well-defined) stochastic process and its 
finite- dimensional laws are given, for any /c G N, by 

(3.9) fi{{xo,...,Xk})= lim P(xf )=xo,...,xf ^=Xfe). 

A' — ^■oo 
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Under fi, X is a (i/, M)-MC, where M is defined as in equation {3.8) and v 
is the invariant probability measure for M and is given by i^(x) = L^Rx, for 

Proof. We begin with (ii) and show that the hmit in (3.9) exists. Let 
/c G N and xq, . . . ,Xk £ E, for any N > k we have 



(3.10) P(xf ) = .0, . . . , = X.) = f n V.,.,. 1 ZlePy 



n 

i=0 



When N goes to infinity, the only significant terms of (V^ '')xk,'xo ^ind 
trace (V^) are those in A^. More precisely, 

(3.11) (V''~%^,xo = Rx,LxX-'+ E flA'A'^-^ 

A'eigenvalue of V^^A 

(3.12) =Rx,LxX-'' + o{X''-') 

as A > I A' I, besides trace(y''^) = A"*^ + o(A^) which leads to 

(3.13) lim F(xf ) = xo, . . . , Xf ) = X.) = n V...... • 

i=0 

Let /i({xo, . . . , Xfe}) = ^""^f"^" 11^=0^ '^x^,Xi+i , we can check that is a proba- 
bility distribution. Indeed, from equations (3.10) and (3.13) 

h h h ^ t--(V^) ^ i^;trace(VA^) 

where the inversion of the sum and the limit is immediate since E is finite. 
We check similarly that the matrix M defined in (3.8) is stochastic. 

Now it is easy to see that the finite-dimensional laws given in (3.9) are 
consistent, the Kolmogorov extension theorem applies and ensures that the 
stochastic process X is well defined. 

Point (i) follows directly from the definition of a cyclic Markov chain. □ 

4. Examples of graphs. We give in this section some examples of results 
that can be obtained by the application of Theorems 2 and 3. In the following 
examples, we only consider oriented lattices with vertices indexed by a subset 
of 1? . The jth line of the graph is the set of vertices with second coordinate 
equal to j. We will see why the restriction of the gas model "to a line" of 
the graph is Markovian. The general approach used is widely inspired by 
the method developed in [4] and [5]. For a given graph G, we first show that 
the assumptions of Theorem 2 are verified for G and the sequence of lattices 
{G^^^)n, that implies that the gas process on converges in distribution 
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to the gas process on G. We then compute the distribution of the gas on a 
Une of G^^^ and interpret it as a cycUc MC by checking that its distribution 
can be written in a multipHcative form as in equation (3.7). Theorem 2 and 
Theorem 3 imply then that the gas process restricted to a Une is Markovian. 
We explain fully the first example and sketch the others. 

4.1. The family of lattices {Cfj)fjc'N. We define in this section a new 
family of lattices. For any finite subset i? of N such that \R\ > 2, we define 
Cr as the lattice with set of vertices indexed by and from each vertex 
(i, j), there are emerging edges from to {i + r,j + 1) for r G R. In 
the following, we always assume that inf (R) = without loss of generality. 
We set R = sup{R). 

Note that J^{o,i} corresponds to the square lattice. If = {0, . . . , n — 1}, 
then Cr = Cn, which corresponds to the family of lattices introduced in [5] 
and detailed in the following. Another example is given in Figure 4. 

Remark 2. For any finite subset R of N, the lattice Cr verifies the 
assumption of Section 2 so the gas model is well defined for any p < p^^f^ 

and since the outdegree of any vertex is equal to \R\, p^^^^^ > l/\R\ > 0. 

For N > n + R, the balls or radius n of C^^^ and of Cr are isomorphic, 
moreover assumption 2 of Theorem 2 holds true with d= \R\, thus the 
finite-dimensional laws of the gas model on Cr converge to the ones on 

■ 

We denote -'^j^'* the iV-tuple that gives the occupation of the gas on the 
jth line of ^C^-* and compute its distribution. The construction of the gas 
model given in Section 2 implies that, when j decreases, (xj^^)jgz is a 
"vertical" Markov chain (with 2^ states) under its stationary distribution. 
Markov chain theory implies that such a distribution is unique (that one of 
the main tool in [5] and [4]). 

For G C [A^], let F(^'^ be the probability that the occupied vertices of a 
line of the graph are exactly those with first coordinate lying in G . In other 




Fig. 4. Example of a DA of size 10 on lattice CR, when R — {0, 1, 4}. 
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words, for a gas occupation distributed according to the gas model 

given in Section 2: 

(4.1) = Pp^'^' [X^^R^ = 1 if and only if i G C). 

Note that the construction of the gas model implies that Pj^^ does not 
depend on a particular choice of j. 

We traduce the fact that iXj^^)j^z,ji is Markovian into recurrence rela- 
tions for . To that purpose, we define for any subset C of N 

M{C)= \J{i + r\reR} 

iec 

and 

Ar{C)= [j{i-r\r£R}, 

where the addition is taken in [N]. Notice that {M{C) x {1}} and {M{C) x 
{ — 1}} correspond respectively to the set of children and of fathers of the 
set {C X {0}} and that \N{C)\ = \M{C)\. We thus obtain the following 
equations: 

(4.2) F^)=f^)"" E 



Following Bousquet- Melon and Conway [5], we check that the probability 
distribution defined by 

\c\ 



^ Zn\1-v) ^ ' 



\XiC)\ 



where 

zn= y: 

CCIN] 

is stationary for equation (4.2). 

To obtain a matrix formulation of that distribution, we consider the ma- 
trix V indexed by the elements of {0, 1}^ and defined by, for a = {sq, . . . , s^_i) 
and r = (ti, . . .,tj{), 

V,^,^ = if (si,...,Sj^_i) / (ti,...,tjj_J; 

and otherwise, 

Vo-,T = ^ ^ — P, if = 9-1^^ there exists r such that Sr = 1 and R — r ^ R, 
1, otherwise. 
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The quantity can then be rewritten as 



N-l 



(4.3) 



with 



Z/7V 



=0 

where 



Cat = (To and 

c7i(A;) = l, iffi + fc-lGC, 




^iv= }__. I I I V„,.„_ 1 =trace(V^^ 

(Tl,...,( 

The expression given for F^'^ in equation (4.3) is in the very same form 
as the statement of Theorem 3. Furthermore it is immediate to check that 
all the coefficients of V^~^ are positive [for p G (0, 1)] which ensures that V 
satisfies the conditions of Theorem 3 by the Perron-Probenius theorem. 

As mentioned in Remark 2, the finite dimensional laws of the gas occu- 
pation on C^^^ converges to those on Cr. We apply Theorem 3 and get: 



Theorem 4. Let X = (^(«, j))(ij)gz2 be the gas process on Cr dis- 
tributed according to Pp^, with p < p^^^^. The stochastic process (5]j)jgN 
defined by Sj = {X{i,0), . . . ,X{i -\- R — 1,0)) is a Markov chain under its 
stationary distribution. 

In other words, the stochastic process {Xi = X{i,0))i^fq is a Markov chain 
with memory R — 1 under its stationary distribution. 

The family of lattices Cn ■ The family of lattices {Cn)n>2 introduced in [5] 
corresponds to the particular case of Cr when i? = {0, 1, . . . , n — 1} (examples 
of DA on £3 and £4 are given in Figure 5). In [5] the GF for DA with one 
source is given as the solution of an algebraic solution of degree at most 
n + 1. 

Theorem 5 (Bousquet-Melou and Conway [5]). The generating func- 
tion Q for DA on Cn with a single source is solution of the following equa- 
tion: 

(4.4) t\l + t)"-^[l + (n + l)g]"+i - [1 + t + (n - l)gY'~\t - 20^) = 0. 

We give some examples of computation obtained by the application of 
Theorem 4 on the lattices £„. In the case n = 2, the computation of the 
eigenvalues and eigenvectors ofV=('^^p ^) constitutes an alternative 



proof of Theorem 3.3 of [10]. 
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Fig. 5. DA on the lattices £3 and £.4 (all the edges are oriented upwards). The cells are 
black and the sources are circled. 



For n = 3, the transition matrix can be given explicitly as (the coefficients 
of the matrix are indexed by the lexicographical order on {0, 1}^) 



/1/A 




(4.5) 



a 




1-1/A 


1 — a 




1 



\ 

p/2X p/2X 


-p/X p/Xj 



where < 



a 



X 



(2-p-A)A' 



l + y/l+4p- V 



For example, we obtain as a consequence of that formula that the generating 
function for DA on £3 with a compact source of size A; > 2 is equal to 



1 - i(Vl-4f-4t2) 



1 - 4t - 4*2 + (1 + 2t) Vl - 4t - 4*2 V 1 + Vl - 4t - 4*2 



-2t 



k-l 



To obtain the GF as the solution of an algebraic equation, we use that in 
[5] , the largest eigenvalue A of V is shown to be solution of 



(4.6) 



n-1 



a"-Ha-i)' 



For n > 4, A cannot be computed explicitly from equation (4.6). Neverthe- 
less, since L and R are eigenvectors associated to A their coordinates can 
be computed in linear time and are polynomial of degree one in A. With 
the condition of renormalization J^^i^i = 1, we obtain that for any free set 
S, the generating function for DA on with source is a rational frac- 
tion and its numerator and denominator are polynomial in A. Moreover, we 
know that A is solution of equation (4.6), which implies that the generating 
function is algebraic in p. 

4.2. The triangular lattice. The triangular lattice, denoted Tri, is defined 
as the oriented graph with set of vertices (i, j) £ Z2 such that i and j have 
the same parity and with set of oriented edges [(«, j), (i — 1, j + 1)], [{i,j), (« + 
l,j + l)] and [{i,j),{i,j + 2)] (see Figure 6). 
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We follow some ideas used in [4] to compute the law of the gas on that 
lattice (note that in [4] the generating function for DA with one source on 
the triangular lattice is obtained by the study of an ad hoc gas model). 
We work on Tri^^^ the cylindric version of the triangular lattice. We keep 
the definition of introduced in equation (4.1), but since a vertex has 
children in the two following lines, we need to define an extension of 
to obtain recurrence relations. Let C and D be two subsets of [N] and X 
a gas model on Tri^^^, we denote F^j^ the probability that the vertices 

occupied in the line (respectively the line 1) of Tri^^-* are exactly the ones 
with first coordinate belonging to C (respectively to D); in other words, for 

C,D€[N], 

j^W _ pTri(^) f ^Tri(^) / • _ i :f ^^u, / £ = and i/2 G C or 
Fc,D -^p [X {^,e)-l if and only if | ^ ^ ^ ^^^^ {^-l)/2GD 

We define M{C) as Uiecl* ~ ^ + 1} which leads to the following recurrence 
relation for C, Z) G [N] such that M{C) nD = 0: 



\c\ 

(4.7) i^gJ=(T^) E i^Ai.(l-p)^-l^(^)^^l 



Notice that since the sum is taken on sets E such that M{D) D E = 0, 
\M{D) U E\ is equal to \J\f{D) \ + \E\. Therefore the distribution given by 

„|CL|D| 

(4.8) Fc,D = ^^Y^ W(C)nD=0 for C,De [N] , 



and where 



C,D 
MiC)nD=0 



is solution to the recurrence relation given in (4.7). 
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Let V = ( ^ Q j , we can rewrite equation (4.8) as 



2N-1 



^""''^ " trace(V^) Q 

(4.9) 

{1, if i is even and i €C, 
1, if i is odd and i ^ D, 
0, otherwise. 

Combining equation (4.9) and Theorem 3 results in the following statement: 



Theorem 6. Let X = (^(i,i))(j,j)eTri be the gas process under Pj", the 
stochastic process 5] = (Sj)jgz defined by 

^ _ ( X{i,0), if i is even, 
* \X[i,l), if i is odd, 

is a Markov chain under its stationary distribution and its transition matrix 
is given by 

(Si = 0|5]o = 0) P(Si = l|5]o = 0)\ /1/A p/A2' 
(Si=0|I]o = l) P(Si = l|So = l)y' ^ 1 



w 



where A = ^+v^^+^p and its stationary distribution is given by 
[P(So = 0), P(5:o = 1)] = [aV(p + X^),p/ip + A^)]. 
Adding up equation (4.8) for all possible D leads to 
(4.10) Fc = ^pl^l(l+p)^-l-^(^)l. 

Setting = (^^^ p) enables equation (4.10) to be rewritten as 



N-l 



(4.11) Fc = TTTivT IT Vxi,xi+i, where = 1 if and only if 2i G C. 

trace (v/^^j 

Again Theorem 3 and equation (4.11) lead to: 

Theorem 7. Let X = (^(j,i))(ij)gTri gas process under Pj''', the 

stochastic process S = (5]j)igz defined by Sj = X{2i,0) is a Markov chain 
under its stationary distribution and its transition matrix is given by 

(Si = 0|5]o = 0) P(Si = l|5]o = 0)~ 
(Si = 0|I]o = l) P(Si = l|5]o = l) 



W 

(4.12) 



1 O^Q O^o 

a, 1 — a. 
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Fig. 7. Examples of sources considered in Proposition on the left and in Proposition 
on the right. 



and its stationary distribution by [P(5]o = 0),P(5]o = 1)] = [-^^f^,-^^^^ 
where 



2p ^ 1 + ^/TT^ 
"°= l + Vl+4p = 2 • 

The Hnk between gas distribution and enumeration of DA given in Propo- 
sition 1 and a simple matrix computation give the following reinterpretation 
in terms of enumeration of DA of Theorems 6 and 7 (see Figure 7 for an 
example of the different sources considered) . 

Proposition 1. (i) Let S = {si, . . . ,Sk} where Si = {xi,€i) be some 
points on the triangular lattice with G {0, 1} and such that di := Xj+i — Xi 
for i £ {1, . . . ,k—l} are non smaller than 2. The GF of DA on the triangular 
lattice with source S is given by 



where a = ■ 

(ii) Xei 5 = {si, . . . , Sfc} where Si = {2xi,0) be some vertices on a line of the 
triangular lattice, such that di := Xj+i — Xi for i £ {1, . . . ,k — 1} are positive 
integers. The GF of DA on the triangular lattice with source S is given by 

^ a, + Oo 1 a, + Oo 

In particular, if 5„ := {(i, 0), i = 1, . . . , n}, we obtain p) = q,^+°q. (1 ~ 

!)"• Then, the GF of DA on the triangular lattice with compact 
sources satisfies 



nil " + ^ + ^p 

That formula was obtained in [8] by combinatorial methods. 
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4.3. The family of lattices Tn. We now study the family of lattices 7^ 
introduced by Corteel, Denise and Gouyou-Beauchamps in [6]. The oriented 
lattice 7^ is a combination of the lattice and the triangular lattice, defined 
as follows: 

• if n = 2/e + 1, the vertices of 7^ are labeled by the elements of Z^. From 
each vertex (i, j) G there are n emerging edges from (i, j) to (i + r, j + 1) 
for —k <r <k and one emerging edge from {i,j) to {i,j + 2); 

• if n = 2k, the vertices are labeled by the elements (i, j) G such that i 
and j have the same parity. From each vertex there are n emerging 
edges from to (i + 2r + 1, y + 1) for —k < r < k — 1 and one emerging 
edge from (i, j) to {i,j + 2). 

The case n = 2 corresponds to the triangular lattice, treated separately 
in Section 4.2 for sake of clarity. In [6], the generating function for DA on 
Tn with a single source is shown to be solution of an algebraic equation 
given explicitly. The proof relies on a combinatorial argument which links 
the generating function for DA on 7^ to that for DA on £„. 

The method used to obtain a stationary distribution for the gas model on 
Tn is very similar to that used in the case of the triangular lattice in Section 
4.2. We keep the same definitions for F^'^ and as those given for the 
triangular lattice and define for C G [N], M{C) as the set, 

• Ujgcl^ + '"i f*^'^ ~^ <r<k} if n = 2A: + 1, 

• Ueci^ + 2r + 1, for -A; < r < - 1} if n = 2k. 

With that new definition of AA(C), equations (4.7), (4.8) and (4.10) still hold 
true. Following the ideas introduced to study £„ given in [5], we define V as 
the square matrix (Vo-,r)(T,T with indices running over {0, 1}"-^ and defined 
as follows. If 0" = (si, . . . ,Sn-i) and r = {t2, ■ ■ ■ ,tn), then 

if (S2, ■ • ■ , Sn-l) / (^2, • • • , tn-l), 

if (s2, ■ • ■ ,Sn-i) = (^2, • ■■,tn-i) and Si = 1, 
if a = r = (0,0,...,0), 
otherwise. 

The stationary distribution of the gas model on a line of 7^^^'' is given by 



(4.13) V, 




(4.14) 



with 



Fn = — n""" V 



where 



o'N = o"o and 

ai{k) = 1, if and only if i + A; — 1 G D, 



^^ = E. nV...,.)=trace(V^). 
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The characteristic polynomial of V, denoted x can be calculated explicitly 
X{x) = (x"" - + 2p) +/ ^ x^\ . 

\ k=0 I 

We rewrite the latter equation as 

X(X) = (p2 - +p2 _ _ (2p + 1))). 

1 — X 

That implies that the dominant eigenvalue A of V satisfies A 7^ 1 and 
(4.15) p2 ^ ^n-i(^ ^^^_X)(\- (2p + 1)). 

We are here in the very same situation as for £„. We can compute explic- 
itly the solutions of equation (4.15) only for n < 4. Nevertheless the same 
arguments as those given for £„, apply and we obtain the following from 
Theorem 3 and equation (4.14). 

Theorem 8. Let X = (X(i, j))(jj)g7-^ he the gas process under FJ". The 
stochastic process (Sj)jgpj defined by 

= {X{2i, 0),X{2{i + 1), 0), . . . , X(2(n - 2), 0)) fori£ N, 

is a Markov chain under its stationary distribution. 

In other words, (X(2i,0))jgN is a Markov chain with memory n—1 under 
its stationary distribution. 
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